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Abstract 


Herein, we construct an explicit modal numerical solver based on the spec- 
tral Petrov—Galerkin method via a specific combination of shifted Cheby- 
shev polynomial basis for handling the nonlinear time-fractional Burger- 
type partial differential equation in the Caputo sense. The process reduces 
the problem to a nonlinear system of algebraic equations. Solving this alge- 
braic equation system will yield the approximate solution’s unknown coef- 
ficients. Many relevant properties of Chebyshev polynomials are reported, 


some connection and linearization formulas are reported and proved, and 
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all elements of the obtained matrices are evaluated neatly. Also, conver- 
gence and error analyses are established. Various illustrative examples 
demonstrate the applicability and accuracy of the proposed method and 
depict the absolute and estimated error figures. Besides, the current ap- 
proach’s high efficiency is proved by comparing it with other techniques in 


the literature. 


AMS subject classifications (2020): 65M60, 11B39, 40A05, 34A08. 


Keywords: Time-fractional Burgers’ equation; Chebyshev polynomials; Petrov— 


Galerkin method; Convergence analysis. 


1 Introduction 


The orthogonal polynomials of the shifted first-kind Chebyshev polynomials 
(SIKCPs) are defined on the interval [—1,1]. Due to their advantageous 
characteristics and value in spectral approaches, they play a vital part in the 
numerical solution of partial differential equations (PDEs). 

Using an orthogonal polynomial-based series expansion, spectral methods 
[11, 25, 19, 20, 10] are numerical approaches that try to approximate the 
solution of a PDE. These techniques are particularly well adapted to the 
orthogonality, clustering, and exponential convergence of the SIKCPs. 

The S1KCPs were denoted as T,,(a%) = cos(n@), where n is the degree of 
the polynomial and @ = cos~1(zx) is the angle between the x-axis and the 
point (2,0) on the unit circle. 

The clustering property of the S1KCPs is one of their main benefits. The 
Chebyshev polynomials are unique among orthogonal polynomial families 
in grouping at the ends of the interval [—1,1]. This clustering property 
is beneficial when approximating functions with boundary layers or steep 
gradients at the endpoints. 

Another feature of spectral methods based on Chebyshev polynomials is 
their exponential convergence. By truncating the series expansion involving 
the Chebyshev polynomials, spectral methods approximate the solution of 


a PDE. As the degree of the truncated series increases, the approximation 
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approaches the real solution exponentially. Chebyshev spectral methods are 


particularly efficient and precise due to their rapid convergence feature. 


Chebyshev spectral methods for numerically solving PDEs require mod- 
eling the solution as a series expansion with S1KCPs as basis functions. By 
projecting the PDE onto the basis functions and solving a sequence of alge- 


braic equations, the series expansion coefficients are determined. 


The number of basis functions (polynomial degree) chosen depends on 
the required accuracy and the issue features. Higher degrees provide more 
precise results but also require more computational resources. In practice, it 


is critical to balance precision and efficiency. 


Chebyshev spectral methods, such as [7, 9, 10], have been effectively used 
for various fractional PDEs, such as elliptic, parabolic, and hyperbolic equa- 
tions. They have proven especially useful in problems with smooth solutions, 


periodic boundary conditions, and unbounded domains. 


Finally, the S1KCPs [27] are valuable tools for numerically solving PDEs. 
They are well-suited for spectral approaches due to their orthogonality, clus- 
tering, and exponential convergence qualities. These polynomials can ob- 
tain accurate and efficient approximations of PDE solutions, making them 
a powerful tool in computational mathematics and engineering. Chebyshev 
polynomials have recently been widely employed to solve several forms of 
differential problems; for example, see [41, 8, 1, 38, 42, 17, 2, 6, 16, 26]. 

We can guarantee that knowledge of the properties and applications of 
S1KCPs in the numerical solution of PDEs is generally established and found 
in various references and textbooks on numerical methods for PDEs. Among 
the significant works in this field are [14, 35, 13]. 


The nonlinear time-fractional equation Burger’s equation is a PDE with 
a fractional derivative in time that incorporates nonlinear convection and 
diffusion factors. It is a variation of the traditional Burger’s equation, a 
simple model for various physical phenomena, such as fluid flow and traffic 
movement. Incorporating fractional derivatives in time allows the system to 


include nonlocal and memory effects. 


The general form of the nonlinear time-fractional Burger’s equation is 


given as [18]: 
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Ou 
aya TOU a hag = 3(,4), O0<a<l, (1) 
u(z,0)=g'(x), O<#<l, 
u(0,t)=9°(t), u(l,t)=9°(t), O<tS1, (2) 


where u(x,t) is the unknown function representing the dependent variable, 
t is time, x is the spatial variable, o is a constant coefficient that controls 
the strength of the convection term, « is a constant coefficient controlling the 


diffusion term, and s(a,t) is the source term. 


Due to the presence of both nonlinear and fractional variables, solving the 
nonlinear time-fractional Burger’s equation analytically is difficult. Numeri- 
cal approaches, on the other hand, can be used to approximate its solutions. 
Many techniques, such as finite difference methods, finite element methods, 
and spectral methods, can be used to deal with fractional derivatives and 


nonlinearity. 


Numerical approaches for solving the nonlinear time-fractional Burger’s 
equation frequently involve griding the spatial domain and using time- 
stepping methods to approximate the temporal derivatives. Fractional dif- 
ference operators and fractional integral transformations can approximate 
fractional derivatives. For further methods that studied other types of frac- 


tional differential equations, see [5, 22, 29, 23]. 


Furthermore, the behavior of the nonlinear time-fractional Burger’s equa- 
tion can exhibit fascinating phenomena, such as the generation of shock 
waves, solitons, and other nonlinear waves. Memory effects are introduced 
by the fractional derivative in time, which might alter the transmission and 


evolution of these nonlinear structures. 


Scientists are constantly researching the properties, analytical solutions, 
and numerical approaches of the nonlinear time-fractional Burger’s equation. 
It is used in a variety of domains, such as fluid dynamics, heat transfer, 
and nonlinear wave phenomena, where the inclusion of fractional derivatives 
in time allows a more precise representation of the system dynamics. The 
interested reader can see the recent diverse numerical methods used to solve 
Burgers’ problem in [15, 31, 44, 36, 24]. 
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The Petrov—Galerkin method [34, 27] is a numerical technique for solving 
PDEs. It is a variation of the more general Galerkin approach, which aims to 
approximate a PDE solution by projecting it onto a finite-dimensional trial 
function subspace. To improve the accuracy and stability of the approxi- 
mation, the Petrov—Galerkin approach incorporates an additional weighting 


function known as the test function or the Petrov function. 


The Petrov—Galerkin approach [43] is especially effective for PDEs with 
specific properties, such as those with convection-dominated terms or singu- 
larities. It solves the usual Galerkin method’s difficulties in capturing precise 


solutions in these challenging settings. 


The Petrov-Galerkin approach [28] has been effectively used for var- 
ious PDEs, including convection-diffusion, Navier-Stokes, and advection- 
dominated situations. It provides a versatile framework for dealing with 


difficult PDEs and is more accurate and stable than the Galerkin technique. 


Overall, the Petrov—Galerkin approach is a strong numerical methodology 
that extends the Galerkin method to meet the difficulties given by specific 
types of PDEs. It improves the accuracy and stability of the approximation 
by introducing Petrov functions as extra weighting functions, making it an 


important tool in the field of numerical PDE solving. 


The structure of this article is as follows: The theory of fractional calculus 
and the relevant properties of Chebyshev polynomials are briefly introduced 
in section 2. A numerical spectral Petrov—Galerkin technique for solving the 
time-fractional Burgers’ type equation is constructed in section 3. Section 4 
discusses the convergence and error analysis of the method. Some numeri- 
cal examples are given in section 5 to illustrate the theoretical conclusions. 


Section 6 contains conclusions. 


2 Preliminaries and essential relations 


2.1 The fractional derivative in the Caputo sense 


Definition 1. [32] The Caputo fractional derivative of order s is defined as 
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1 


Diu(x) = yea i (a —y)™ sul (y)dy, s>0, 2>0, (3) 


where m—l<s<m, meEN. 


The following properties are satisfied by the operator D} for m—1< 
s<m, me€EN, 


D*c=0, (cis a constant) (4) 


x 


0, if me No and m< [sl], 
—— (5) 
rete, if mENo and m>{s], 


where N = {1,2,3,...$ , No = {0} UN and the notation [a] denotes the 


ceiling function. 


2.2 An account on the SIKCPs 


Let T(x) be the SIKCPs defined in the interval [0, 1] by T7 (x) = Tj(2 2-1). 
These polynomials can be defined as [8, 38] 


ey eR (H 1) * 2?* (m+ k 1)! 
ZT; aS (m— bh) eg 0, (6) 


satisfying the following orthogonality relation with respect to the weight func- 


tion w(x) = CE [8, 38]: 
1 
F, iba) 2 (a) T* (2) de Fin Sov ns (7) 
0 
where 
a, if m=0, 
hm = (8) 
Ss, if m>0, 
and 
1, ifm=n, 
0, ifmFn. 


The recurrence relation of T*,(z) is 
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Tin¢i(@) = 2 (2@—1) Ty, (@) — Tri (2), (10) 


where Tj (a) = 1 and Ty (#2) =2a—-1. 


Moreover, the inversion formula is [8, 38] 


Jj — 91- 23 ( )! * : 
wl = 2'~*7 (23) a jap ee) 720, (11) 
p= at 
where 
4, ifm =0, 
i (12) 
1, otherwise. 


The following relation between T;‘ (x) and U;*(x) is correct 
DT} (x) = 2iUj;_,(x), for alli > 1. (13) 
The following linearization formula is valid 


1 
Tea) = 5 (Uz,,.() +UZ_,(z)), for all r,s > 0. (14) 


Corollary 1. [1] For every positive integer q, the gth derivative of T7 (x) can 


be expressed in terms of their original polynomials as 


& 
Q 


D'T; (x) = dj.p,q Tp (2), (15) 
0 
) 


p= 
(j+p+q) even 


where 


a 52? Peg (9) ng —p-4) 
dina = EG apap + Mw . 


and €, is defined in (12). 


Lemma 1. [4] For all nonnegative integers m and n, the following lineariza- 
tion formula holds for the S1KCPs: 


Te, (a) T2(2) = 5 (Thal) + Ton—ni(2)) - (17) 


Lemma 2. Let 7 and 7 be any two nonnegative integers. The moments’ 


formula for the SIKCPs are given by 
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ne P(2j +1) ' 
la a Dd, GG+]— s+ it jeeey ve 


s=i— 


Proof. Multiplying both sides of (11) with T;*(x) and direct use of Lemma 


1, we get the desired result. 


Remark 1. The following relation is satisfied: 


1 
[ PU @) To) dr = oi0, (19) 
0 
where U;‘(x) is the shifted Chebyshev polynomials of the second kind and 


m, if (¢—r)even,i>r, 
Our = (20) 
0, otherwise. 


Remark 2. The following relation is satisfied 


1 
' x w(x) Te (x) T3 (x) dx = kis, (21) 
0 
where 
if |s — i] = 2, 
if |s — i] =1, 


ifs—i=0,i>1, 


(22) 
ifs=7=0, 


ifs=i=1, 


, otherwise. 


3 Petrov—Galerkin approach for the time-fractional 


Burgers’ equation 


In this section, we consider the following time-fractional Burgers’ equation 
[33]: 


O° x(a, t) 


Ox(@,t) 4, x(a.) 
v 
ote 


+ x(z,t) an 22 


=S(a,t), O<a<l, (23) 


subject to the following initial and boundary conditions: 
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xy(2,0) = g(a), O<a <1, 
x(0,)=(t), x0,Hj=G(t), O<t<1, (24) 


where W is the kinematic viscosity and S(x,t) is the source term. 


Now, to proceed with our proposed Petrov—Galerkin approach, we will 


use the following transformation: 
x(2,t) = u(x,t) + T(2,t), (25) 
where 
T(z, t) = (1-2) (x(0, 4) — x(0,0)) + & (x(1, 4) — x(1,0)) + x(z,0), (26) 


to convert (23) governed by the conditions (24) into the following modified 


equation: 
O° u(a, t) Ou(a,t) OY (a,t) 
Rye + Ula, t) —37— + ula, t) a 
2 
+ T(a,t) MD _ yg FU) _ pi), 0<a <1, (27) 


governed by the following homogeneous conditions: 


u(z,0)=0, O<a<1, 


u(0,t) =u(1,t)=0, O0<t<1, (28) 
where 
7 8° T(z, t) 8 T(z, t) oT (at) 


Therefore, instead of solving (23) governed by (24), we can solve the modified 


equation (27) governed by the homogeneous conditions (28). 


3.1 Trial functions 


Consider the following basis functions: 


Xj (a) = Tio(a) — T(x), 
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oj (t) = tT; (t). (30) 


Remark 3. The polynomials A(x) can be written alternatively in the fol- 
lowing form: 


Ai (a) = —8a (1 — 2) UF (2). (31) 


The orthogonality relations of \F(x) and ¢j(t) are given by 


: * * 1 a 
[rn @) MO ape = BF ba (32) 


and 


, * * 1 = 
i nD OO ge Bee = hm bn (33) 


where h,, is defined in (8). 


Lemma 3. For all nonnegative integers m and n, the following linearization 


formula, holds: 


«(gy #8) __ sy Lose ae 
(0) SE) = sp pal) + 24 +2) Ut mal) — 6+ 2) UE ng) 
+ 4UG m—1(@) — (24+ 2) UE maa (@) + (6+ 2) Ui mia(2). 

(34) 


Proof. We express D F(x) as a combination of U(x) via (13). Then we 


linearize A* (x) D X* (a) using (14), we get the desired result. 


Lemma 4. For all nonnegative integers m and n, the following linearization 
formula holds: 
* * t * * 


Proof. The proof of this lemma is a direct result of Lemma 1. 


Theorem 1. The first and second derivatives of W*,(a) can be expressed 


explicitly as 


dM(2) SG. 
ae = Do si T; (x), 
j=0 
@P¥(@2) Goo 
ant = 27 iy (x), (36) 
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where 


j+1, ifi+1=3, 


2, if («-— 7) odd,j > 0, 
oi =4 

is if (i — j) odd,j =0, 

0, otherwise, 


8(3i(i+2)—j2 +4), if G—j)even,i>j, 
75,4 = 95 
0, 


otherwise, 


1, iff =0, 
65 = (37) 


$, otherwise. 


Proof. Relations (36) can be deduced after using Corollary 1 when g = 1,2 
along with \*(x) defined in (30), then collecting like terms and rearranging 


the summations. 


3.2 Petrov—Galerkin solution for the time-fractional 


Burgers’ equation 


Now, one may set 


0% (Q) = span{A# (x) 6; (t) : i,j =0,1,...,.N}, 
AN (Q) = {u € OX (Q) : u(x, 0) = u(0,t) = u(1, t) = 0}, (38) 


where 2 = [0,1]?. Then any function uw (x,t) € AX (Q) may be written as 


N WN 
uN (x,t) = > > ca AT (@) HO. (39) 


i=0 j=0 
The application of Petrov-Galerkin technique [39] is used to find uw (x,t) € 
AN such that 


(2 uN (a, t) 


ur (x 
Fe ATO) + (Wen mane) 


w(a,t) 
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OL (et) ‘ 
+ ( w (a, £) ———"—, Ti (a) TS (t) 
( ( ) Joc 
Ou (a, t) x * 
ee a ae tet Ty (x t 
oO? uN (a, t) * * 
wU Tr(a) Te (t 
( Ox? (2) Ts ) t) 


= (F(o,t),T(2) TH))uiey OS SSN 


(40) 


where T;* (a) T; (t) is the test function and w(x,t) = w(x) w(t). 


Therefore, (40) can be written after using the definition of w (x,t) (39) 
as 


N N N WN 


N WN 
s s Cij Gir bj, as s S s s Cmn Cij dm. ir ere 
1=0 j=0 m=0 n=0 i=0 7=0 

N N N WN N N N WN 


a S- S- S- S- Cmn aij dm. ir hn j,s + S- > > x bmn Cig dm, ir An. j,s 


m=0 n=0 i=0 7=0 m=0 n=0 i=0 j=0 


N oN 
Vy Soa; Ziv Qj,s = Frs, OSr, sc N, (41) 
i=0 j=0 


where a;; and b;; are determined from the following relations: 


 - if ae fe.) ee 1 1 
“=a | P Bae A () 65 (8) (cd —a)e? B 7c ee 
"81h; Jo Jo iA cae cL. ce ae 

(42) 
Also, 
; : d F(t) 
tir = OCA) TED ate te (SE re) 


dX* 
dinar = (Aila) 7 (2) so Ie (HOGO. TO a 


a? \3(x) mais 
Lie = ( x) ) Ty (2) ie ’ Qj,s =. (65 (t), i ()) sq ) 
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and Fis = (f(a, t), fie (x) a (t) (wt) : 


Theorem 2. The elements gj, bjs, dm,ir, hn,j,s, Zi,r, and Q;,5 are given 


by 
-2, if r-i=0, 
Gir = a if r—1=2,, 
0, otherwise, 
3 2, =I 9; 3 a 
bj. =m(-1T (3a) 4F3 1], 
5, s—a+2,s—a+2 
Amir =-1 Oi—-m-3,r 1 (2% oF 2) Oi—m-1,r — (i T 2) Oi—-m+1,r +1 Oitm—l1,r 


_ (24 + 2) Oitm-+1,r + (i + 2) Oi+m43,r> 


1 
hn,j,s = 5 (Kj4n,s ae Kj—n,s) ’ 


4m (3i(i+2)—1r? +4), if @—r)even,i>r, 


Zi => 
0, otherwise, 
Ber? ifs—j=1, 
=, ifj—-s=1, 
Qj, = = ; ; (43) 
ai ifs—j=0, 
0, otherwise. 


Proof. The elements g;, can be obtained after using the orthogonality re- 
lation (7) and the definition of A¥(a) defined in (30), then collecting like 


terms. 


To find the elements b;,,, based on the definition of Caputo’s fractional 


derivative (3) and the power form of T(x) (6), one can write 


I 52? 4 1)M—DI-FG4R—D! flows tne 
2 aby Hlca eka | Te (tt 8 w(t) dt 
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25 vee 1)J-*D(k + 2)0 (7 + kD (k — a + 3) 


. (44 
a P(2k + ITC k+1)[(kK-s—at+2)(k+s—a+2) Ga 
The last relation can be summed to give the following result: 
3 2, Id; 3 a 
bjs =m (-1/P TP (5 = a) 4f3 hs (45) 


$,-8—-a+2,s—at2 


The elements dy,;,- can be obtained after using Lemma 3 along with Remark 
il: 

Similarly, the elements h,;,, can be obtained after using Lemma 4 along 
with Remark 2. 

The elements Z;,, can be obtained after using Theorem 1 along with the 
orthogonality relation (7), collecting like terms, and rearranging the summa- 
tions. 


Finally, the elements Q;,, can be obtained after using Lemma 2 along 


with the orthogonality relation (7) and doing some computations. 


Remark 4. The inner product of (u(x,t), v(2,t))o¢@,1) > (W(), U(@)) aca), and 
(u(t), v(t)) sq) are defined as 


(u(x,t), 0(2,4)) ce.) = | | lett dude 
(u(2), 0(2)) ace) = [ uae) v() @(@) de, 
(u(t), v(t) ac) = | u(t) v(t) w(t) dt. (46) 


Remark 5. The (N + 1) x (N + 1) nonlinear system of equations in (41) 
in the unknown expansion coefficient cj; can be solved through a suitable 


numerical solver such as Newton’s iterative technique. 


4 Error bound 


Herein, we give an upper bound of the absolute errors (AEs) using Lagrange 


interpolation polynomials. 
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Algorithm 1: Coding algorithm for the proposed technique 

Input ¥W, g(x), Ci(t), Ga(t), a, and S(z, t). 

Step 1. Using transformation (25) to convert the nonlinear TFCE 
(23)—(24), into modified equation (27)—(28) 

Step 2. Assume an approximate solution wv (z, t) as in (39). 

Step 3. Apply Petrov-Galerkin method to obtain the system in (41). 

Step 4. Use Theorem 2 to get the elements of 9; , bjs, dm,r hn,j,s, Zi,rs 
and Qj,s- 

Step 5. Use FindRoot command with initial guess 
{c;; = 10~*4, 1,7: 0,1,..., N}, to solve the system (41) to get cj. 

Output wu (z,t). 


Let uw (x,t) € AN (Q) be the best approximation of u(z,t); then, the def- 


inition of the best approximation enables us to write the following inequality: 


u(x,t) — uN (2, t)|loo < |lu(a,t) — uv (x, t)\|o0, for all v% (a, t) € AN (9). 
(47) 

Moreover, the previous inequality is also true if @ denotes the interpolating 
polynomial for u(x,t) at points (x;,t;), where x; are the roots of A*(x), while 


t; are the roots of $7 (t). 


Using similar steps as in [12, 40], one has 


NHu(nt) Fy ONHu(e,w) Ty 
= N = : 2 rm 
u(a, t) v (x, t) OaNtl (N +1)! I @i) + NGI (N +1)! He tj) 
PNP ui) 7 7 
daN+1 AtNt1 ((N + 1)!)?2 1G x) [[@-+). 
i=0 j=0 
(48) 
where 1,4), 1, fl € (0, 
Now, 
alia! II To (a = @4)II 
Ber < ee 
||u(x,t) — v" (2, t)|loo S fe neti wth : a Ty 
se anme [Ort wa 4)| [TT Lo(t= tlle 
(a,t)EQ a (N +1)! 
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Pn ult B) 


OxrNt+1 OtNtl 


(rtjen 
NTI o(@ — #i)lloo I TTLo(t = tlle 
(V+)? : 


(49) 


Since u is a smooth function on 2, then there exist three constants (1, 2, 
and €3 such that 


ON+! u(ar, t) ON+1 u(a, us) 
ax_|—~—~ | < 4, ax bo, 
(ater OgNr1 aoe (@.t)e0 OtNt1 alia 
PN*? u(h, ft) 
max |-——~—— | < és. 50 
(xtjen|OaN+1 OtNt1 | — eu) 
To minimize the factor || To (e — 2X;)|loo, let us use the one-to-one mapping 


a = 4 (z+1) between the intervals [—1, 1] and [0, 1] to deduce that 


N N 
: 1 
wi€(0.1] eae Il = i) 7 ee 1] te | ra Il a C _ i) 
1\ N+ N 
= (5) ge he 


N+1 
1 r 
= = 1 
(5) Zi ay 1] AL. mieT , ? ) 


where \y = 2% is the leading coefficient of Ay—1(z) = Ty41(z) — Ty_-1(2) 


and z; are the roots of Ay_1(2). 


Also, the factor || Ik 3G —t;) ||. may be minimized by using the one-to- 
one mapping t = 3 (¢ + 1) between the intervals [—1, 1] and [0,1] to deduce 
that 


(i) 
On 


N { N+1 
min max [[@-+t) = (5) min se ; (52) 


t; €[0,1] te[0,1] 7x0 €; €[-1,1] te[- 


where dn = 2"? is the leading coefficient of dv (é) = a Ty (t) and ¢; 
are the roots of @y(t). 


Since 


max, [Aw—a(2)| = [Zn+1(1)| + [Zy—1(0)] = 2 
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ace lon+i(#)] = lon4i(1)| = 1, (53) 


te[-1, 


then, inequalities (50), (51), (52), and (53) enable us to get the following 


desired result: 


|lu(x, t)—w (x, t)|loo < £1 = (5)% ; Gr ie G 2N41 
; ee XN (N+UL bw (N+U)E 7 Aw bn (N “oF 
54 


which represents an upper bound of the AE. 


5 Illustrative examples 


Example 1. [83, 37] Consider the time-fractional Burgers’ equation of the 


em a ult Ou(ast) —. Pula) 
“~ula,t u(az,t u(x,t) 
ofa + u(z, t) az 2 Ax? = g(z,t), (55) 
subject to the following initial and boundary conditions: 
u(z,0)=0, O<a<1, 
u(0,t) =u(1,t)=0, O<t<1, (56) 


where g(z,t) = Toa) ?-“sin(nx)+7t4 sin(m x) cos(r x) +277 t? sin(r 2), 
and u(x,t) = t? sin(7 2) is the exact solution of this problem. 

Table 1 gives a comparison of AE between our method and the method 
in [33] at a = 0.7 and a = 0.8. Figure 1 shows the maximum absolute error 
(MAE) at different values of N when a = 0.5. Table 2 gives a comparison 
of L-error between our method and methods in [37, 33] at different values 
of a. Figure 2 shows the AE (left) and approximate solution (right) at 
a=0.9,N =14. 


Example 2. [33] Consider the time-fractional Burgers’ equation of the form 


O° u(a, t) 
ote 


anaes oul 2 ue Dandy, (57) 


subject to the following initial and boundary conditions: 


u(z,0)=0, O<a<1, 
u(0,t)=2?, u(l,t)=—#, O<t<1, (58) 
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Table 1: Comparison of AE for Example 1 


a=0.7 a=0.8 


x Method in [33] Our method at M=14 Method in [33] Our method at N = 14 


0.1 2.18701 x 10-5 5.55112 x 10-17 8.10773 x 107-7 5.55112 x 107-1” 
0.2 4.20505 x 107-5 1.11022 x 10~16 1.57153 x 107 1.11022 x 10~16 
0.3. 5.88630 x 10-5 2.22045 x 10-16 2.22713 x 1076 1.11022 x 10—16 
0.4 7.06958 x 10-5 3.33067 x 10-16 2.71560 x 10-6 
0.5 7.61346 x 10-5 2.22045 x 10-16 2.97247 x 10~® 
0.6 7.41772 x 10-5 3.33067 x 10~16 2.94204 x 10~® 
0.7 6.44989 x 10-° 4.44089 x 10716 2.59370 x 10-6 1.11022 x 10-16 
0.8 4.76900 x 10~° 3.33067 x 10716 1.93833 x 10-6 1.11022 x 10-16 
0.9 2.53573 x 10-5 5.55112 x 107!” 1.03759 x 10-® 5.55112 x 107-1” 
10% 
2 1078 
10-8 
4 6 8 10 12 14 
N 


Figure 1: The MAE of Example 1 at a = 0.5. 


Table 2: Comparison of Loo-error of Example 1 


a Method in [37] at N = 27, M=2' Method in [33] at N= 27, M=2'? Our method at N = 14 


0.2 5.47402 x 107% 4.76511 x 107% 3.33067 x 10716 
0.3 5.46367 x 1075 4.75486 x 107° 3.33067 x 10716 
0.4 5.44879 x 107% 4.74015 x 10~® 2.22045 x 10-16 
0.5 5.42862 x 10-° 4.72029 x 107° 7.77156 x 10716 


2 


where g(x, t) = pagtay t?% cos(w x) — wt* sin(m x) cos(m x) + 1° t? 


cos(7 x), 
and u(x,t) = t? cos(7 2) is the exact solution of this problem. 

Table 3 reports the AE at N = 12 and N = 14 when a = 0.5. Figure 
3 shows the MAE at different values of N when a = 0.9. Table 4 gives a 


comparison of L,.-error between our method and method in [33] at different 
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Figure 2: The AE (left) and approximate solution (right) for Example 1 at a = 0.9, N 


14. 


Approxi 


imate solution 
a 


values of a. Figure 4 shows the exact and approximate solutions at a 


190 


0.4,N = 14. 
Table 3: The AE of Example 2 at a = 0.5 
N=12 N=14 
x t=0. t= 05 (a0 t=0.1 t=0.5 t=0.9 
0.1 5.05238 x 10-16 1.26565 x 107 4.10366 x 107'4_ 6.50521 x 107! 6.93889 x 10718 4.16334 x 1071” 
0.2 1.22949 x 10-' = 3.07462 x 10-!4 9.96425 x 10-4 2.60209 x 10-18)»: 1.88778 x 107!" 5.55112 x 10-17 
0.3 1.43115x 107! 3.57214x 10-4 1.14991 x 10-18_—-8.67362 x 107! 2.77556 x 10-17 1.11022 x 10-16 
0.4 1.04387 x 10-18 2.59792 x 10- 8.32667 x 10-4 1.95156 x 10-'8 1.73472 x 10-17_—-1.38778 x 10-17 
0.5 1.88161 x 107! 1.16453 x 107!” 2.50372 x 10-17 2.69252 x 107! 1.10588 x 10717 5.19085 x 1071” 
0.6 1.04431 x 107! = 2.60071 x 107 8.32112 x 10-4. 2.81893 x 1078 ~—- 3.81639 x 107-17 1.24901 x 10716 
0.7 1.43071 x 10719 3.56937 x 107 1.15047 x 10718 8.67362 x 10718 1.38778 x 10717 0 
0.8 1.22949 x 10-1 3.07601 x 10-4 9.95871 x 10-4 2.60209 x 10-18 — 2.77556 x 10-17 0 
0.9 5.05238 x 10-16 = 1.26427 x 10- 4.10644 x 10-4 6.50521 x 107! 2.08167 x 10-17_—-9.71445 x 10-17 
A076 
10° 
uw 
< 
= 
40-1" 
49° 


Figure 3: The MAE of Example 2 at a = 0.9. 


14 
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Table 4: Comparison of Loo-error of Example 2 


Method in [33] 


a N=2?,M=2" M=2",N=2" Our method at N = 14 


0.2 1.06427 x 1075 9.66513 x 1076 4.44089 x 10-16 
0.3 1.06304 x 10-5 9.60898 x 10-6 4.44089 x 10-16 
0.4 1.06126 x 10-° 9.53096 x 107° 3.33067 x 10716 
0.5 1.05886 x 1075 9.42186 x 1076 3.33067 x 10716 


Exact solution pi ale solution 


Figure 4: The exact and approximate solutions of Example 2 at a = 0.4,N = 14. 


Example 3. [33] Consider the time-fractional Burgers’ equation of the form 


O° u(z, t) Ou(z,t) 6? u(z,t) 
ae t - t 59 
Fae tala t) A — REA = o(2,), (59) 
subject to the following initial and boundary conditions: 
u(z,0)=0, O<a#<1, 
u(0,t)=2#?, u(lt)=et?, O<t<1, (60) 


where g(a,t) = Ta) ef t27-% 4 ¢4 ¢2® _ #? e? and u(x,t) = t? e” is the exact 
solution of this problem. 

Table 5 gives a comparison of AE between our method and method in 
[33] at a = 0.9. Table 6 gives a comparison of L..-error between our method 
and methods in [33] at different values of a. Figure 5 shows the AE (left) 
and approximate solution (right) at a = 0.3, N = 12. 
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Table 5: Comparison of AE for Example 3 at a = 0.9 


Method in [33] 


z« N=200,At=0.0005 N =80,At=0.001 Our method at N = 12 


0.1 6.18496 x 10-9 2.23568 x 107-7 5.55112 x 10717 
0.2 7.14787 x 107° 3.99993 x 10-7 1.11022 x 10716 
0.3 3.49621 x 107° 5.30327 x 1077 2.22045 x 10716 
0.4 3.99695 x 10-9 6.14693 x 1077 2.22045 x 10716 
0.5 1.42437 x 1078 6.52199 x 10-7 2.22045 x 10716 
0.6 2.56063 x 10-8 6.40791 x 107-7 0 

0.7 3.55155 x 107-8 5.76962 x 10-7 2.22045 x 10716 
0.8 3.98360 x 1078 4.55240 x 1077 

0.9 3.17808 x 10-8 2.67253 x 10-7 


Table 6: Comparison of L.o-error of Example 3 


Method in [33] 


a N=2',M=2'!' M=27,N=24 Our method at N = 12 


0.2 5.69465 x 1077 1.97043 x 10~° 2.22045 x 107-16 


0.3 5.66088 x 1077 1.99857 x 107° 2.22045 x 107-16 


0.4 5.61387 x 1077 2.03716 x 107° 2.22045 x 107-16 


2.22045 x 10-16 


2.09231 x 10° 


0.5 5.54967 x 10-7 


Approximate solution 


Figure 5: The AE (left) and approximate solution (right) for Example 3 at a = 0.3, N = 
12. 


Example 4. Consider the following time-fractional Burgers’ equation of two 


dimensional 
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age + tw) ( ie °° oy Ox 
FU e.yyt) 
7 g(x,y, t), (61) 


subject to the following initial and boundary conditions: 


u(z,y,0)=0, O<a,y<1l, 
u(0,y,t)=t?e 4, u(l,y,t)=t?e 4%, O<y,t<1, 


i 00=fe, dnljsre, 0<e0< 1, (62) 
where g(z,t) = ray" — 2t?e7-¥ and u(z,t) = t? e*-¥ is the exact 


solution of this problem. 


Table 7 illustrates the AE at different values of t when N = 6 and a = 0.2. 
Figure 6 shows the AE at different values of t at N = 6 and a = 0.2. 


Table 7: The AE of Example 4 at a = 0.2 


r=y t=0.3 t=0.6 t=0.9 
0.2 5.50355 x 1074 2.19933 x 10-18 4.94896 x 10718 
0.4 4.41578x 107! = 1.7732 x 10718 4.01347 x 10-8 
0.6 4.40518 x 107! 1.75651 x 107'8 3.92865 x 1078 
0.8 5.50272x 10-4 2.19789x 10-8 4.94157 x 10-18 


Error At t=0.1 Error At t=0.5 
el, ee, fl , 


a —— wf yee: 


Figure 6: The AE for Example 4 at a = 0.2, N = 6. 
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6 Concluding remarks 


In this study, using a special combination of shifted Chebyshev polynomial 
bases, we built an explicit modal numerical solution based on the spectral 
Petrov—Galerkin technique to handle the nonlinear time-fractional Burger- 
type PDE in the Caputo sense. The procedure reduces the issue to a set of 
nonlinear algebraic equations. Numerous important Chebyshev polynomial 
characteristics were reported, along with various connection and linearization 
equations that were mentioned and verified. All components of the resultant 
matrices were also elegantly assessed. Additionally, studies of convergence 
and error were created. The applicability and accuracy of the suggested 
technique were illustrated through a number of illustrative instances, which 
also show the absolute and anticipated error rates. As an expected future 
work, we aim to employ the developed theoretical results in this paper along 
with suitable spectral methods to treat some other problems, for instance, 
(21, 3, 30]. All codes were written and debugged by Mathematica 11 on an HP 
Z420 Workstation, Processor: Intel (R) Xeon(R) CPU E5-1620 - 3.6 GHz, 16 
GB RAM DDR3, and 512 GB storage. 
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